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ABSTRACT 

We propose the indirect search for the CP violating phase in the long baseline 
oscillation experiment, in which two scenarios of the neutrino mass hierarchy 
are discussed in the three family model. The CP violatig phase effect is clearly found 
in the scenario: the LSND data plus the atmospheric neutrino deficit. The phase 
dependence of the oscillation probability is explicitly shown by using typical parameters 
of the K2K experiment. The matter effect is negligibly small. In order to select the 
scenario of the neutrino mass hierarchy, the measurement of the '^r oscillation is 
also proposed in the long baseline experiment. 
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1 Introduction 



Neutrino flavor oscillations provide information of the fundamental property of neu- 
trinos such as masses, flavor mixings and CP violating phase. In these years, there is 
growing experimental evidences of the neutrino oscillations. The exciting one is the 
atmospheric neutrino deficit |T|~0 as well as the solar neutrino deficit ||^. Super- 
Kamiokande also presented the large neutrino flavor oscillation in atmospheric neu- 
trinos. Furthermore, a new stage is represented by the long baseline(LBL) neutrino 
oscillation experiments [0. The flrst LBL reactor experiment CHOOZ has already re- 
ported a bound of the neutrino oscillation 0, which gives a strong constraint of the 
flavor mixing pattern. The LBL accelerator experiment K2K is planned to begin 
taking data in the next year, whereas the MINOS and ICARUS [jlO| experiments 
will start in the flrst year of the next century. Those LBL experiments will clarify 
masses, flavor mixings and CP violation of neutrinos. 

Some authors have already discussed possibilities of observing CP violation 
in the LBL experiments by measuring the differences of the transition probabilities 



between CP-conjugate channels ||12||, which originates from the phase of the neutrino 
mixing matrix, such as — >■ Vr and — > z7^. However, the direct measurement is 
very difficult in the planned LBL experiments since the magnitude of its difference is 
usually expected at most 0.01 and the difference of energy distributions of neutrino 
beams and z7^ disturbs this measurement in the order of (9(0.01). 

In this paper, we propose the indirect search for the CP violating phase in the 
LBL experiments combined with results of the short baseline(SBL) experiments. The 
SBL experiments give important constraints for neutrino mixing angles. The tentative 
indication has been already given by the LSND experiment ll^, which is an accelerator 



experiment for z/^ — > u^.i^p^ Ve)- The CHORUS and NOMAD experiments have 
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reported the new bound for Vr oscillation, which has already improved the E531 



result [|T^. The KARMEN experiment [jT6| is also searching for the t'^ — > t'e(^/x ^ 



oscillation as well as LSND. The Bugey ^\ and Krasnoyarsk [T^ reactor experiments 



and CDHS [jl9[ and CCFR ||2^ accelerator experiments have already given bounds for 



the neutrino mixing parameters as well as E776 ||21 |. 

Taking account of those data, we study the CP violatig phase effect of the z/g 
oscillation in the LBL experiments. Bilenky, Giunti and Grimus have already 
discussed the bound of this transition probability, which is independent of the CP 
violatig phase. By investigating the phase effect of the bound, we provide a possibility 
to observe the CP violation in the neutrino oscillation indirectly. 

2 Mass and mixing patterns and LBL experiments 

Let us start with discussing the recent results of atmospheric neutrinos at Super- 
Kamiokande [Q, which suggests — * Vj. oscillation with the large mixing. Since the 
CHOOZ result excludes the large neutrino oscillation of — > z/g, the large mixing 
between and Vr is a reasonable interpretation for the atmospheric z/^ deficit. Our 
starting point as to neutrino mixings is the large — Vr oscillation with 

Am^,^ - 5 X 10-=^eV2 , sin^ 2^atm ^ 1 , (1) 

which constrain neutrino oscillations in the LBL experiments. The considerable large 
oscillations are expected for the z/^ process in K2K 0, MINOS and ICARUS 
lOl experiments. This large oscillation hardly depends on other mixings and the CP 



violating phase. On the other hand, the z/^ v^. oscillation is small and depends on 
the phase and other mixing angles, which are constrained by the SBL experiments and 
the LBL reactor experiments. 

In order to give the formulation of the neutrino oscillation probability for the LBL 



i^fi experiment, we should discuss at first the mass and mixing pattern of neu- 

trinos. Other possible evidences in favor of neutrino oscillations are the solar neutrino 
and the accelerator experiment LSND. Three different scales of mass-squared differ- 
ences cannot be reconciled with the three family model of neutrinos unless one of data 
is disregarded. Therefore, one should two possible scenario|2^: (1) "sacrifice solar 
neutrino" and (2) "sacrifice LSND". By stretching the data, the other scenario is 



still available as discussed by Cardall and Fuller This scenario appears to emerge 
naturally as the most likely solution to all oscillation evidences, however, contradicts 
with the zenith angle dependence of multi-GeV atmospheric data in the Kamiokande 
and Super-Kamiokande experiments 0. So, we do not discuss this scenario in our 
paper. If one would like to sacrifice no data, the sterile neutrino should be included in 
the analyses, which is out of scope in our paper. 

In the scenario (1) "sacrifice solar neutrino", we take Amg-^ ~ IS.vn^^ ^ Am^gND 
and Am|2 — Am^^.jjj(m3 ^ ^ vti}). The natural pattern of the neutrino mixing 



matrix U 25 



IS 



f/ - ha c s , (2) 
\e4 -s c) 

where c = cos^^ ~ and s = sin 9 ~ are fixed by taking account of eq.(|l|), 

and small values e^'s are constrained from the SBL experiments. We do not consider 
the mass hierarchy ^ m2 ^ rrii, which is unfavourable in constructing the neutrino 
mass matrix without fine-tuning in the case of the large mixing P3 . 



In the scenario (2) "sacrifice LSND", we take Am^g — ^f^li — Am^^^^ and 
Am^i — Ami, (solar neutrino). For the MSW small angle solution |2^, the mixing 



pattern is the same one in eq.(|2|). Another pattern is equivalent to an exchange of the 
first and second columns in the mixing matrix. In the case of the just-so solution and 
the MSW large angle solution, the mixing matrix pattern is somewhat complicated. 
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Taking account of the data of the CHOOZ experiment, Bilenky and Giunti [^] have 
obtained a typical mixing pattern, on which we will comment later. 

Now, in terms of the standard parametrization of the mixing matrix U 

/ C13C12 C13S12 Si3e"*'^\ 

U = — C23S12 — S23Si3Ci2e"^ C23C12 — S23Si3Si2e*'^ S23C13 , (3) 
\ ■S23S12 — C23Si3Ci2e*'^ —-523012 — £23^135126*'^ C23C13 / 

where = sin6'jj and c^- = cos 6ij are vacuum mixings, and is the CP violating 
phase, the neutrino oscillation probability for the LBL — > experiment is given by 

P(z/^ Ve) ~ 2\Uei\^\U^i\^ - 4c23S23Si2Si3COS(/)sin^ — -j^ 2C23S23S12S13 siu Scp 

(4) 

with 

\Uel\^\U^a\^ ~ C23S12 + S23S13 + 2c23S23Sl2Sl3COS(/) , (5) 

for the scenario (1), and 

P{.Vf, i^e) ^ 4S23S13 sin^ 2C23S23S12S13 sin Sep (6) 

for the scenario (2), where 

Am?r,L Am'^nL Aml.L , , 

S^cp = sin ^ + sin ^ + sin (7) 

2E 2E 2E ^ ' 

The oscillatory term Sep is approximately reduced to sin(Am23i^/2i?) for the scenario 
(1), because sin(Am^2-^/2-£') and s\Ti{/S.ml^L/2E) can be replaced by the average value 
in the LBL experiments. For the scenario (2), Sqp is suppressed such as Sep ~ 0.01, 
so the second term in eq.(^ can be neglected. 

Thus, the probability in eq.(|^) strongly depends on the CP violating phase, while 
the one in eq.(||) hardly depends on it. The CP violating phase effect could be found 
in the LBL — > z/g experiment if the scenario (1) is true. On the other hand, if the 
scenario (2) is true, it is impossible to determine the CP violating phase by z/g. 



3 Constraints from reactor and accelerator experi- 
ments 

In our analyses, S23 — l/-\/2 is fixed, but S12, S13 and (p still remain as free patrameters. 
Those mixing parematers are constrained from other experiments. For the scenario 
(1), the SBL experiments give strong constraints because the distance L corresponds 
to the neutrino mass scale (9(leV). Let us consider constraints from disappearance 



experiments. The Bugey |jT^ and Krasnoyarsk |18| reactor experiments and CDHS |19| 
and CCFR accelerator experiments give bounds for the neutrino mixing parameters 
at the fixed value of Am|^ [^. In the mixing pattern of eq.(0), the constraint for Uei 
is given by using 



P(Z7, ^ 17,) ~ 1 - 4|?7eir(l - It/eiH siu^ 



4E ' 

where ^m\^ ^ Arrigg is used. Thus, one obtains the bound |f/eiP > o-e due to non- 
observation of neutrino oscillations. For example, we get a, = 0.984 for Arrig^ = 2eV^. 
This condition is expressed in terms of mixing parameters as follows: 



■^12 "I" -^la '^12'^13 — -^e — 1 O'e ■ (9) 

Constraints from the appearance experiments are given by LSND, E776 and CHO- 
RUS/NOMAD. If the excess of the electron events at LSND are due to the neutrino 
oscillation, the —>■ Vf. oscillation transition probability is equal to 0.31 ±0.09 ±0.05% 
[rB|. The plot of the LSND favored region of Am^ vs sin^26' allows the mass region 



to be larger than Am^ ~ 0.04 eV^. However, Bugey|17| has already excluded a mass 
region lower than Am^ = 0.27 eV^. On the other hand, the constraint by E776 |21| 
has excluded a region larger than Am^ = 2.3 eV^. Thus, the mass squared difference 
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Am|^ is obtained in the range PO]: 



Am^i ~ 0.27 eV^ ~ 2.3 eV^ . (10) 

The recent CHORUS experiment for has given the bound sin^ 26'chorus < 

1.3 X 10^'^[|l^. Those data lead to constraints of mixing angles by using following 
equations: 

P(^.,-^.e)^4|[/,,n[/,ipsin2^^, 

P{u, v^) ^ ^\UM^A^ sin' . (11) 

Using eq.(II), we get the condition for the LSND data 

C23"Si2 + "523"Si3 + 2C23S23S12S13 COS = = ^ siu^ 26'lsnd , (12) 

where Le is given experimentally for fixed Arrig^. 

For the CHORUS bound of Vt-, we also obtain a condition 

■S23'5l2 + C23-5?3 " 2C23S23S12S13 COS < = ^ siu^ 26'CHORUS , (13) 

which is a weak bound at present since sin^ 26'chorus > 0.1 in the case of Amg^ < 2eV^. 
Actually, mixing angles which satisfy eqs.(^) and (|I2|) always satisfy this bound. 

The first LBL reactor experiment CHOOZ [0 also constrains mixing angles by the 
following formula: 

P{V^ -.V^)..!- 2\U^,\\l - If/.iH - A\UMUe,?sm' • (M) 

Taking account of S12 ^ 1 and S13 <^ 1, the probability reduces to 

P(I7,-. 77,) ^1-2(4 + 4) . (15) 

Since CHOOZ reported sin^2^cHooz < 0.18 for large Am^, we get a constraint + 
-Sis < 0.18/4, which is a loose bound compared with the one in eq.(|^). 
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Next we consider the scenario (2), in which the SBL experiments do not constrain 
mixing angles because the highest mass difference scale is (9(0.01eV^). The constraint 
follows only from the CHOOZ experiment [0, in which the probability is expressed as 

P(F. -17.) ^ l-4|f/.3p(l-|f/e3p)sin^^^ ^ l-4s?3(l-4)sin^^^ . (16) 

By using the CHOOZ result sin^ 26'chooz < 0.12 for Am^ ~ 5 x lO^^eV^, we get only 
a constraint s^g < 0.03. The MSW small angle solution of the solar neutrino gives 
another constraint: Si2 = 0.02 ~ 0.06. 

4 Phase dependence of P{i^^ I'e) 

In the scenario (1), eq.(|D is written as 

P(z/^-z/e) ~2|[/einf/^i|'-4c23S23Si2Si3sin(^^jsin(^^-0j , (17) 

where the first term in the right hand side is fixed by the LSND data. It is remarked 
that C23S23S12S13 determines this oscillation probability if is fixed. So, the maximal 
value of C23S23S12S13 gives the maximal or minimal P{iyfi — > i^e), which depends on 0. 

We show how to extract the maximal value of S12S13 from eqs.(§) and ([T^) , taking 
■523 = 1/ a/2 for simplicity. Those conditions turn to 

4 + 4-44 <5e , (18) 
"S12 + "S13 + 2Sl2Sl3COS0 = 2Le . (19) 

The maximal value of S12S13 under the condition in eq.(0) is obtained by an easy 



algebraic manipulation as follows: 



S12S13 < —^—-L, . (20) 
1 + cos m 



On the other hand, we get an equation as to S12S13 by a subtraction of eq.(|18D and 
eq.dH: 

44 + 2Si2Si3 COS + (5e - 2Le) > . (21) 
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Taking a positive solution of S12S13, we get 

■Sl2'Sl3 < — cos — \J COS^ 4> — {Be — 2Le) . (22) 

The true maximum of S12S13 should satisfy both conditions of eqs.(pOD and (^21) . In 
conclusion, we obtain the maximal S12S13 as follows: 

S12S13 < — ^ -Le for COS0 > -1 + ^(1 + Jl- Be) , 

1 + COS (p Be ^ 

S12S13 < -COS0- ^COS2 0- (5e -2Le) for COS < -1 + ^(1 + \/ 1 - Be) . (23) 

By using this maximal value, we can predict the allowed region of the oscillation prob- 
ability of eq.(p!7D for the fixed value of (f). We show the allowed region of Piv^i —>■ Ve) 
versus in fig.l, where Am|2 = 5 x 10"'^ eV^ and typical parameters of the K2K ex- 
periment L = 250Km, E = 1.3GeV are taken. We also take Amg^^ = 2eV^, which 
leads to Be = 0.0183 and Le = 6.5 x 10^'^. The allowed region is shown between the 
solid curve and the horizontal dashed line. It is noticed that this result is consistent 
with the one in ref. [^. Since the condition of eg. (p3D is a general one, it is useful even 
if the LSND data is replaced by the KARMEN |16| data in the future. 

As seen in fig.l, the probability increases steeply in the region > 100° and to the 
maximal value around 0.02. Thus, the observation of -P(z^^ Ve) larger than 0.005 
indicates the large CP violating phase although the CP conserved limit = 180° is 
still allowed. In the case of = 180°, the observed Piv^ Ve) fixes both Si2 and S13. 
Then, one can predict the CP conjugated process Ve for arbitrary neutrino beam 

energies. Thus, the measurement of the absolute value of Piv^ Ve) is a crucial test 
of the CP violation. In fig.l, we have also shown the oscillation probability of 17^ Ve 
by taking same values of parameters L and E. 

In the scenario (2), substituting a bound 5^3 < 0.03 into eq.(||), we obtain 

P{^t. ^ i^e) < 0.06 , (24) 
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which is much larger than the maximal value 0.02 in the scenario (1). In this case, the 
CP phase dependence is not found. 

We have taken the MSW small angle solution of the solar neutrino in the scenario 
(2). In the case of the just-so solution and the MSW large angle solution, the mixing 



matrix pattern is somewhat complicated. However, the mixing matrix given in ref.pT 
hardly changes our result with the MSW small angle solution. 

The observation of 0.02 < P{Vfj, —>■ Ve) < 0.06 means that the scenario (2) is true, 
that is to say, the mass difference scale Am^ = (9(leV^) is rejected. This is the same 
issue in ref.[0. What can one say if its observed probability is lower than 0.02? The 
measurement of the Ve oscillation can select a scenario because the CP violating 
phase effect is found in the scenario (1). An alternative method is to measure the 
z/g — > oscillation in the LBL experiment. For the scenario (1), the probability is 
different from the —>■ v^. oscillation as follows: 

P(z/e-^Z/.) ^2|f/einf/,i|2 + 4C23S23Sl2Sl3Sin(^^jsm(^^ + 0j , (25) 

while the probability is the same one for the scenario (2) in the case of S23 = l/"\/2- 
Thus, the v^. — * z/^ oscillation can select a scenario of the mass hierarchy. 

It may be useful to remark that the observation of P{yp, —>■ Vg) larger than 0.06 
indicates new physics. 

5 Matter effect 

Although the distance travelled by neutrinos is less than 1000 Km in the LBL experi- 
ments, those data include the background matter effect which is not CP invariant p2| . 
In particular, the matter effect should be carefully analyzed for the scenario (1) since 
the numerical results strongly depend on the CP violating phase. 
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Hamiltonian in matter H{x) in weak basis is 

/ml/2E 
H = U\ ml/2E \W + \ \ , (26) 
V ml/2E^ 

where 

a = y^GpUe , (27) 

with a constant matter density p = 2.34 g/cm^. For antineutrinos, effective Hamiho- 
nian is given by replacing a —a and U ^ U* . 

Minakata and Nunokawa have shown the perturbative treatment of the matter 



effect which is rehable under the mass scale we consider. In the scenario (1) with 



rria ~ 1112 ^ mi, the transition probability is given: 

P{y^. ^ ^e) = 2\Vei{a)\^\V^M? - 4Re[Ve2(a)V;*3(a)V;*2(a)V;3(a)] sin^ I/23 

+ 2Im[V;2(a)V;3(a)V;*2(«)^,.3(a)]sinJ23 , (28) 

where 

K.(a)=f/..-E2ai5;^%^, h, = -^^ + m.?~\U.,?)aL. (29) 

The oscillation probability turns to a complicated equation, but it is easy to estimate 
the matter oscillation effect by using the parameters E = 0{lGeV), a = O{10~^^eV) 
and L = 250Km. We have found the leading matter corrections for the vacuum os- 
cillation to be dominated by terms of 2ai?/Am|i, 2aii^sf2/Am|2, 2aEs1^/Am'^2 
aL{sl2 — s1^)- Those corrections are at most 0(10"'^) for the vacuum oscillation prob- 
ability of z/^ — > Ue in the case of Am|i ~ leV^. This result is consistent with the 



previous numerical one [0. Thus, the matter effect is negligible small. 

6 Conclusions 



We have proposed the indirect search for the CP violating phase in the LBL exper- 
iments combined with results of the SBL experiments in two scenarios of the mass 
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hierarchy. Taking account of those data, we have investigated the CP violatig phase 
effect of the — > oscillation in the LBL experiments for the scenario (1): the LSND 
data plus the atmospheric neutrino deficit. For the scenario (2): the atmospheric neu- 
trino deficit plus the solar neutrino deficit, the phase dependence is negligibly small. 
As an example, we have shown the phase dependence of the oscillation probability by 
using typical parameters of the K2K experiment for the scenario (1). The matter effect 
has been found to be negligibly small. This phase effect on the probability provides a 
possible method to observe the CP violation in the neutrino oscillation indirectly. In 
order to select the scenario (1) or (2), the measurement of the — > Vt oscillation is 
useful in the LBL experiment as well as V^^Vf.. Even if our predicted magnitudes are 
out of sensitivity in the K2K experiment, we can expect the experiments of MINOS 
and ICARUS. 



12 



References 



[1] K. S. Hirata et al, Phys. Lett. B205 (1988) 416; B280 (1992) 146; 
D. Casper et al., Phys. Rev. Lett. 66 (1991) 2561; 
R. Becker-Szendy et al, Phys. Rev. D46 (1992) 3720. 

[2] NUSEX Collaboration, M. Aglietta et al., Europhys. Lett. 15 (1991) 559; 

S0UDAN2 Collaboration, W. W. M. Allison et al., Phys. Lett. B391 (1997) 491; 

Frejus Collaboration, K. Daum, Z. Phys. C66 (1995) 417; 

MACRO Collaboration, S. Ahlen et al., Phys. Lett. B357 (1995) 481. 

[3] Y. Fukuda et al., Phys. Lett. B335 (1994) 237. 

[4] GALLEX Collaboration, Phys. Lett. B327 (1994) 377; B388 (1996) 384; 
SAGE Collaboration, J. N. Abdurashitov et al., Phys. Lett. B328 (1994) 234; 
Homestake Collaboration, Nucl. Phys. B38 (Proc. Suppl.) (1995) 47; 
Kamiokande Collaboration, Nucl. Phys. B38 (Proc. Suppl.) (1995) 55; 
Super-Kamiokande Collaboration, Y. Fukuda et al, [hep-ex/9805021| (1998). 



[5] Super-Kamiokande Collaboration, Y. Fukuda et al, [hep-ex/9805006| (1998); 



T. Kajita, invited talk at Neutrino 98, A - 9 June, 1998, Takayama, Japan. 
[6] S. Parke, |hep-ph/930427T . 



[7] The CHOOZ Collaboration, M. Apollonio et al, Phys. Lett. B420 (1998) 397. 

[8] Y. Suzuki, Proc. of Neutrino 96, Helsinki, June 1996, edited by K. Enqvist et al., 
p. 237 (World Scientific, Singapore, 1997). 

[9] S. G. Wojcicki, Proc. of Neutrino 96, Helsinki, June 1996, edited by K. Enqvist et 
al., p. 231 (World Scientific, Singapore, 1997). 



13 



[10] ICARUS Collaboration, P. Cennini et al., LNGS-94/99-I, May 1994; 

F. Pietropaola, invited talk at Neutrino 98, A - 9 June, 1998, Takayama, Japan. 

[11] M. Tanimoto, Phys. Rev. D55 (1997) 322; Prog. Theor. Phys. 97 (1997) 901; 
J. Arafune and J. Sato, Phys. Rev. D55 (1997) 1653; 
J. Arafune, M. Koike and J. Sato, Phys. Rev. D56 (1997) 3093; 
M. Koike and J. Sato, |hep-ph/9707203| ; 

H. Minakata and H. Nunokawa, Phys. Lett. B413 (1997) 369; Phys. Rev. D57 
(1998) 4403 ; 

S. M. Bilenky, C. Giunti and W. Grimus, |hep-ph/9712537 . 



[12] S. Pakvasa, in Proc. of the XX International Conference of High Energy Physics, 
Madison, Wisconsin, USA, 1980, ed. L. Durand and L. G. Pondrom (AIP, New 
York 1980) Part 2, p.ll65; 

V. Barger, K. Whisnant and R. J. N. Phillips, Phys. Rev. Lett. 45 (1980) 2084; 
Phys. Rev. D23 (1981) 2773; 

S. M. Bilenky and F. Nidermaler, Sov. J. Nucl. Phys. 34 (1981) 606. 

[13] LSND Collaboration, C. Athanassopoulos et al., Phys. Rev. Lett. 75 (1995) 2650; 
77(1996)3082; Phys. Rev. C54(1996) 2685; |nucl-ex/9706006| (1997); 
J. E. Hill, Phys. Rev. Lett. 75 (1995) 2654. 

[14] NOMAD Collaboration, J. Altegoer et al., CERN-EP-98-057 (1998); 
CHORUS Collaboration, E. Eskut et al., Phys. Lett. B424 (1998) 202; 
O. Sato, invited talk at Neutrino 98, A - 9 June, 1998, Takayama, Japan; 
J.J. Gomez-Cadenas, invited talk at Neutrino 98, A - 9 June, 1998, Takayama, 
Japan. 

[15] E531 Collaboration, N. Ushida et al., Phys. Rev. Lett. 57 (1986) 2897. 

14 



[16] KARMEN Collaboration, Nucl. Phys. B38 (Proc. Suppl.) (1995) 235; 
B. Zeitnitz, invited talk at Neutrino 55, 4 - 9 June, 1998, Takayama, Japan. 

[17] B. Achkar et al., Nucl. Phys. B434 (1995) 503. 

[18] G.S. Vidyakin et al., Pis'ma Zh.Eksp. Thor. Fiz. 59 (1994) 364; 
JETP Lett. 59 (1994) 390. 

[19] CDHS Collaboration, F. Dydak et al, Phys. Lett. B134 (1984) 281. 

[20] CCFR Collaboration, I. E. Stockdale et al., Phys. Rev. Lett. 52 (1984) 1384; 
Z. Phys. C27 (1985) 53; 

CCFR Collaboration, K. S. McFarland et al., Phys. Rev. Lett. 75 (1995) 3993. 

[21] E776 Collaboration, L. Borodovsky et al., Phys. Rev. Lett. 68 (1992) 274. 

[22] S. M. Bilenky, C. Giunti and W. Grimus, Phys. Rev. D57 (1998) 1920. 

[23] A. Yu. Smirnov, Proc. of 28th International Conference on High-energy Physics 
(ICHEP 96), Warsaw, July 1996, edited by Z. Ajduk and A.K. Wroblewski, Vol.1, 
p.288 (World Scientific, Singapore, 1997). 

[24] C.Y. Cardall and CM. Fuller, Phys. Rev. D53 (1996) 4421. 



[25] M. Matsuda and M. Tanimoto, |hep-ph/97lT294 . 



[26] S.P. Mikheyev and A.Yu. Smirnov, Yad. Fiz. 42 (1985) 1441; 
L. Wolfenstein, Phys. Rev. D17 (1987) 2369; 

E.W. Kolb, M.S. Turner and T.P. Walker, Phys. Lett. B175 (1986) 478; 

S.P. Rosen and J.M. Gelb, Phys. Rev. D34 (1986) 969; 

J.N. Bahcall and H.A. Bethe, Phys. Rev. Lett. 65 (1990) 2233; 



15 



N. Hata and P. Langacker, Phys. Rev. D50 (1994) 632; 

P.I. Krastev and A. Yu. Smirnov, Phys. Lett. B338 (1994) 282. 

[27] S. M. Bilenky and C. Giunti, ^ep-ph/980220T . 

[28] Particle Data Group, R. M. Barnett et al., Phys. Rev. D54 (1996) 94. 

[29] S. M. Bilenky, A. Bottino, C. Giunti and C. W. Kim, Phys. Lett. B356 (1995) 
273; 

G. L. Fogh, E. Lisi and G. Scioscia, Phys. Rev. D52 (1995) 5334. 

[30] N. Okada and O. Yasuda, Int. J. Mod. Phys. A12 (1997) 3669; 

S. M. Bilenky, C. Giunti and W. Grimus, Eur. Phys. J. CI (1998) 247; 

G. L. Fogh, E. Lisi and G. Scioscia, Phys. Rev. D56 (1997) 3081. 

[31] H. Minakata and Y. Shimada, Astropart. Phys. 8 (1998) 193. 

[32] T. K. Kuo and J. Pantaleone, Phys. Lett. B198 (1987) 406; Phys. Rev. D35 
(1987) 3434; 

P. I. Krastev and S. T. Petcov, Phys. Lett. B205 (1988) 84; 

H. W. Zaglauer and K. H. Schwarzer, Z. Phys. C40 (1988) 273. 



16 




Figure 1: 

Fig. 1: Predicted bound of the — > Ve oscillation in the K2K experiment for the 
scenario (1), which is denoted by the thick sohd curve. The CP conjugated I?^ — > 
one is shown by the thin solid curve. The horizontal dashed line denotes 2|t/eip|?7^ip, 
which is fixed by the LSND data. Here, Amli = 2eV^ and Am|2 = 5 x 10"^eV^ are 
taken. 
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